Anisotropic Hubbard model on a triangular lattice — spin dynamics in HoMnOs 
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The recent neutron-scattering data for spin-wave dispersion in HoMn03 are well described by 
an anisotropic Hubbard model on a triangular lattice with a planar (XY) spin anisotropy. Best fit 
indicates that magnetic excitations in HoMnOs correspond to the strong-coupling limit U /t >~ 15, 
with planar exchange energy J = 4t 2 /U ~ 2.5meV and planar anisotropy AU ~ 0.35meV. 



There has been renewed interest in correlated elec- 
tron systems on triangular lattices, as evidenced by re- 
cent studies of antiferromagnetism, superconductivity 
and metal- insulator transition in the organic systems k — 
(BEDT - TTF^X^ the discovery of superconductivity 
in Na x Co02.yH 2 0^ the observation of low-temperature 
insulating phases in some VS-adlayer structures such 
as K on Si[lll]^ and quasi two-dimensional 120° spin 
ordering and spin-wave excitations in RbFc(Mo04)2 
(Refs. 5,6) and the multifcrroic materials YMnC>3 and 
HoMn0 3 A£ 

Recent neutron-scattering studies of the multifcrroic 
material H0M11O3 have revealed a non-collinear 120° an- 
tifcrromagnctic (AF) ordering below Tn ~ 72 K of the 
5 = 2 Mn 3+ spins arranged in offset layers of two- 
dimensional (2D) triangular latticed Measurements of 
the spin wave dispersion were found to be well described 
by a nearest-neighbour Heisenberg AF with exchange en- 
ergy J = 2.44 meV and a planar anisotropy D = 0.38 
meV at 20 K. No discernible dispersion was observed in 
the out-of-plane direction, indicating primarily 2D spin 
dynamics. 

Recently spin- wave excitations in the 120° AF state of 
the Hubbard model on a triangular lattice were studied 
within the random phase approximation (RPA) in the 
full U ranged The spin wave energy in the large U limit 
was shown to asymptotically approach the correspond- 
ing result for the Quantum Heisenberg antifcrromagnct 
(QHAF), thus providing a continuous interpolation be- 
tween weak and strong coupling limits. However, com- 
peting interactions and frustration were found to signifi- 
cantly modify the dispersion at finite U, resulting in van- 
ishing spin stiffness at U ~ 6 and a magnetic instability 
at U ~ 7 corresponding to vanishing spin-wave energy 
at wave vector q M = (7r/3,7r/v3). The sharp fall-off of 
com near U ~ 7 provides a sensitive indicator of fmite-C/ 
effects in the AF state. Indeed, recent high-resolution 
neutron-scattering studies of the spin- wave dispersion in 
the square-lattice S=l/2 AF La2Cu04 have revealed no- 
ticeable spin-wave dispersion along the MBZ edge^ as- 
sociated with finite-?/ double-occupancy effectsiii 

In this brief report we extend the spin-wave analy- 
sis to include planar spin anisotropy and show that the 
neutron-scattering data for HoMnOs are also well de- 
scribed by a Hubbard model on a triangular lattice, thus 
providing a microscopic description of the most essen- 
tial features of spin dynamics in the 120° AF state of 
HoMnC>3, including the spin- wave dispersion and energy 



scale. We examine the behaviour of spin-wave anisotropy 
gap with spin anisotropy and also suggest a sensitive mea- 
sure of finite-C/, double occupancy effects. 

However, the Mn spin planar anisotropy is treated here 
only at a phenomenological level, equivalent to the ef- 
fective anisotropy DSf z included in recent investigations 
using spin models£& In a detailed study of single-ion 
anisotropy and crystal-field effects in the layered rare- 
earth cuprates R2CUO4 (R=Nd,Pr,Sm), the magnetic 
behaviour (including spin-reorientation transitions) has 
been attributed to coupling of Cu with the rare-earth 
magnetic subsystem which exhibits a large single-ion 
anisotropy resulting in preferential ordering of rare-earth 
moments along specific lattice directions^ It has been 
suggested that the anisotropy of Mn spins, its observed 
temperature dependence, and the reorientation transi- 
tions in HoMnOs also originate from a similar anisotropic 
exchange coupling with the rare-earth Holmium^ result- 
ing in magnetic behaviour as seen in layered rare-earth 
cuprates, where frustrated interlayer coupling allows for 
weaker, higher-order interactions to control the magnetic 
structure. Especially relevant for non-collinear ordering, 
the anisotropic exchange (Dzyaloshinski-Moriya) inter- 
action D.(S; x Sj) originating from spin-orbit coupling 
has been suggested as responsible for the clamping of 
ferroelectric and antiferromagnetic order parameters in 
YMn0 3 ii£ 

Hund's rule coupling responsible for the 5 = 2 spin 
state of Mn +++ ions and crystal-field splitting have 
also not been realistically incorporated here. However, 
these realistic details do not qualitatively affect the 
spin-rotation symmetry and spin dynamics, as discussed 
below. Hund's rule coupling in the generalized Hub- 
bard model considered here is maximal as inter-orbital 
Coulomb interaction for parallel spins is dropped com- 
pletely. Including an inter-orbital density-density inter- 
action Vq and an intra-atomic exchange interaction Fq 
favouring parallel-spin alignment (Hund's rule coupling), 
the more realistic orbital Hubbard modeli^i^ 

i,7<7',£7,<7' 

(1) 

i,7<7 / ,fT^<7 / 

remains spin-rotationally invariant under a global rota- 



2 



tion of the fermion spin S^ 7 = ^L^^^ (where = 

( a i'y^ a tyi.) ^ S ^ ne f errmon held operator), even if orbitals 
7 are identified with the Mn orbitals (t2 g ,e g ) resulting 
from crystal-field splitting of the atomic 3d orbitals. As 
the intra-atomic exchange interaction Fq is much larger 



than the spin excitation energy scale 



within a 



U+F ' 

strong-coupling expansion) , all fermion spins on a site are 
effectively coupled, yielding a composite quantum spin S 
and a corresponding multiplication by factor 2S to ob- 
tain the effective spin-wave energy scale. Therefore, or- 
bital multiplicity does not change the Goldstone-mode 
structure, and the spin-dynamics energy scale in the or- 
bital Hubbard model is essentially determined by t and 
U cS = U + F . 

As a simplest extension to phenomenologically include 
spin-space anisotropy, while retaining only the relevant 
energy scales t and J7 e ff, we consider the generalized Af- 
orbital Hubbard modeli& 



^2 \ "> f t 



(2) 



on a triangular lattice with nearest-neighbour (NN) hop- 
ping between sites i and i+S. Here 7, 7' refer to the (ficti- 
tious) degenerate Af orbitals per site. The factor i is in- 
cluded to render the energy density finite in the Af — ► 00 
limit. The two correlation terms involve density-density 
and exchange interactions with respect to the orbital in- 
dices. The Hartree-Fock (HF) approximation and Ran- 
dom Phase Approximation (RPA) are of 0(1) whereas 
quantum fluctuation effects appear at higher order within 
the inverse-degeneracy expansion and thus 1 /Af , in anal- 
ogy with l/S for quantum spin systems, plays the role of 

h. 

The key feature of spin-rotation symmetry of the gen- 
eralized Hubbard model is highlighted by writing the two 
interaction terms as 



int 



Af^ 



U 2 -Ut 



Af 



in terms of the total spin operator 



(3) 



(4) 



where if>\ = (a\ 



L.i). An Ising (uniaxial) anisotropy 



is obtained for U\ > U2, a planar (XY) anisotropy for 
U2 > Ui, and full spin-rotation symmetry for U\ = U2- 

As appropriate for H0M11O3, we consider the case U2 > 
U\ corresponding to preferential ordering of spins in the 
x — y plane in spin space and an anisotropy gap for out- 
of-planc excitations^ Magnetic excitations were analyzed 
in terms of a Heisenberg model with a similar anisotropy 
term in Ref. [7]. At the HF level, the interaction term 



for orbital 7 then reduces to 



(5) 



where the sclf-consistently determined mean field Aj = 
E^2(Si 7 ')HF lies in the x — y plane in spin space. The 
HF sublattice magnetization depends only on [7 2 and is 
determined from the self-consistency condition (Sq,)hf = 
z(kj ^ItI^j I) a in terms of the HF states |k, I) on sub- 
lattice a, exactly as for the isotropic case£ In the strong 
coupling limit the sublattice magnetization (S a ) « 1/2 
at the HF level, and is reduced by about 50% (in the 
isotropic case) due to quantum fluctuations, as found in 
different calculations cited in Ref. [9]. 

We consider the 120° ordered AF state on the trian- 
gular lattice, and examine transverse spin fluctuations 
about the broken-symmetry state. At the RPA level, the 
magnon propagator reduces to a sum of all bubble dia- 
grams where the interaction vertices involving Sf x and 
Sf y appear with interaction U2, whereas those involving 
Sf z with interaction U\. Introducing a planar spin rota- 
tion, so that spins are oriented along the x' direction for 
all three sublattices, we obtain for the transverse spin- 
fluctuation propagator 



[x(q, 



J a/3 



IVW)] 



2[C/][x°(q,u,)] 



(6) 



in the 2<g>3 spin-sublattice basis of the two transverse spin 
directions [i, v = y',z' and the three sublattices a, (3 = 
A. B. C . The sublattice-diagonal interaction matrix 



[U] 



U 2 1 








(7) 



in the y' , z' basis, and the bare particle-hole propagator 

[x°(q^)G 



4 ^ 

k.t.m 



E k-d,m ~ E k.l + U E k,m - E k-q,l - w 



(8) 



involves integrating out the fermions in the broken- 
symmetry state. In the particle-hole matrix elements 



q, m 



c M |k, l) c 



(9) 



of the rotated spins, the spin orientation angles (f> a in the 
fermion states |k, I) are transformed out. The numerical 
evaluation of [x°(q, ^)] in terms of the HF-level AF-state 
energies and amplitudes is exactly as for the isotropic 
case studied earlier^ 

The spin-wave energies w q are then obtained from the 
poles 1 — Aq(wq) = of Eq. (6) in terms of the eigenvalues 
A q (o;) of the matrix 2[U] [x°(q, <*>)]. As w q corresponds 
to spin 1/2, it is scaled by the factor 2S for arbitrary 
spin SH As expected for planar anisotropy, there is only 
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FIG. 1: (color online) Spin-wave dispersion for the three 
modes calculated from Eq. (5) with J = 4t 2 jJJ ~ 2.5meV and 
AU w 0.35meV, along with neutron scattering data points for 
HoMnOs at 20K from Ref. [7]. 



one Goldstone mode corresponding to planar rotation of 
spins, and the two out-of-plane modes become massive 
with an anisotropy gap w gap . 

Spin-wave dispersion w q for the three modes is shown 
in Fig. 1 along two symmetry directions in the magnetic 
Brillouin zone (MBZ), along with the neutron-scattering 
data for HoMn03 at 20 K from Ref. [7]. The anisotropic 
Hubbard model provides a remarkably good description 
of the spin dynamics. We hnd that best fits with the 
magnetic excitations in HoMnOs are only obtained in 
the strong-coupling limit U/t >^ 15, with a planar ex- 
change energy J = 4t 2 /U ~ 2.5meV and anisotropy 
U-2 — Ui = AU ~ 0.35meV, the individual values of t 
and U not being resolvable within experimental resolu- 
tion. To estimate the order of magnitude of the ratio U/t, 
if we nominally take U = lcV, we obtain t = 25meV and 
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FIG. 2: (color online) The ratio wi/w2,3 at wave vector qju 
provides a sensitive indicator of finite-C/, double-occupancy 
effects, shown for different values of spin anisotropy AU. 
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FIG. 3: Spin- wave anisotropy gap corresponding to out-of- 
plane fluctuations shows a ~ (AU/U) 1 / 2 behaviour with spin 
anisotropy. Here U/t = 16. 



U/t = 40. Spin-wave dispersion calculated in the inter- 
mediate coupling regime cannot be fitted with the neu- 
tron scattering data, indicating no evidence of finite-C/ 
double occupancy effects, as discussed below. 

For the square-lattice S — 1/2 AF La2Cu04, finite- 
U, double-occupancy effects associated with higher-order 



4 



(t 4 /U 3 ) spin couplings are manifested in noticeable spin- 
wave dispersion along the MBZ boundarypifliii Similarly, 
for the isotropic triangular-lattice AF, the ratio u>i/uj2.3 
of the non-degenerate and degenerate spin- wave energies 
at wave vector = (n/3,Tr/y/3) is a sensitive mea- 
sure of the U/t ratio^ which asymptotically approaches 
2/VlO = 0.632 in the strong-coupling (U/t — > oo) 
limit, decreases monotonically with U/t, and eventually 
vanishes at U/t m 7~ Variation of oji/uj2,3 with U/t 
is shown in Fig. 2 for different values of anisotropy 
AU/t. Neutron-scattering data for HoMn03 shows that 
w i/ w 2,3 ~ llmcV/16meV w 0.7, which yields a lower 
bound (~ 15) on the ratio U/t from Fig. 2. 

The spin-wave anisotropy gap u> gap /t, corresponding to 
out-of-planc fluctuations with q = 0, varies as (AU /U) 1 ^ 2 
with spin anisotropy (Fig. 3), which translates to w gap oc 
y/JAU, the geometric mean of the two energy scales. 

In conclusion, the anisotropic Hubbard model provides 
a simple extension to phcnomcnologically include spin- 
space anisotropy and study spin excitations in the full 
range from weak to strong coupling. Comparison of spin- 



wave dispersion with RPA calculations allows for a quan- 
titative determination of the effective anisotropy, in ad- 
dition to highlighting any finite-?/, double occupancy ef- 
fects as seen in La2Cu04. The recent neutron-scattering 
data for spin-wave dispersion in HoMn03 are found to 
be well described by an anisotropic Hubbard model on 
a triangular lattice with a planar (XY) spin anisotropy. 
We find that the twin constraints w gap /wi.2 ~ 1/3 and 
w i/^2,3 ~ 2/3 in the neutron scattering data cannot be 
satisfied in the intermediate coupling regime, and best 
fit indicates that magnetic excitations in HoMn03 cor- 
respond to the strong-coupling limit U/t >~ 15, with 
J ~ 2.5meV and AU ~ 0.35meV. The ratio U\/w2,z of 
the non-degenerate and degenerate spin-wave energies at 
wave vector qj\/ = (tt/3, n/y/S) is suggested as a sen- 
sitive measure of finite-?/, double occupancy effects in 
a triangular-lattice antifcrromagnct. Finally, in view of 
the formal resemblance with the orbital Hubbard model 
(1), our RPA calculation provides a step forward towards 
investigating magnetic excitations using realistic models 
including Hund's rule coupling, crystal-field splitting etc. 
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